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Abstract 

We consider the problem of small data global existence for a class of semilinear wave 
equations with null condition on a Lorentzian background (R 3+1 , g) with a time dependent 
metric g coinciding with Minkowski metric outside the cylinder { (t, x) \ \x\ < R}. We show 
that the small data global existence result can be reduced to two integrated local energy 
estimates and demonstrate these estimates in the particular case when g is merely C 1 close 
to the Minkowski metric. One of the novel aspects of this work is that it applies to equations 
on backgrounds which do not settle to any particular stationary metric. 

1 Introduction 

In this paper, we study the behavior of solutions to the Cauchy problem 

= 7=B 9 « (9 aP V=Gd P 4>) = F^, deb), G = det(g), 
[0(0, x) = €(j) O {x),d t (l)(0,x) = ecj)i(x) 



(1) 



on a Lorentzian manifold (R 3+1 ,g) with initial data 4>o(x), <fii(x) <G C^°(R™). The nonlinearity 
F is assumed to satisfy the null condition with respect to the Minkowski metric 

F(0,0) = 0, dF(0,0) =0, 

F(0, 90) = A^d a 4>d^ + O(|0| 3 + |50| 3 ), (2) 
in which A alB are constants and A al3 ^ a ^/s = whenever £ 2 = £ 2 + £f + £§. 

The corresponding Cauchy problem in Minkowski space 



□0 = F{<t>,d<j>, <9 2 0), 

0(0, a;) = e0 o (a;),9 t 0(O,x) = e0i(x) 



(3) 



has been studied extensively. A classical result is the short time existence for large initial data [7] , 
[29] . The problem of long time behavior of the solutions with small data has been investigated 
by Fritz John [8], [10]. In [8] , application of the standard energy methods led to a lower bound 
on the time existence T > O(-), while in [9], it was shown that any nontrivial C 3 solution of 
the equation 

□ = 2 

with compactly supported initial data blows up in hnitc time. For general nonlinearity, an 
almost global existence with a lower bound on the time existence T > Cexpj^} was established 
in [11]. These results and the corresponding results in higher dimensions [13], [17], [26] rely 
on the decay properties of solutions to a linear wave equation on R™ +1 . 

A remarkable progress in understanding the problem of small data global existence in three 
and higher dimensions has been achieved by S. Klainerman [15]. His approach relied on the 
vector field method, which connects the symmetries of the flat M n+1 with the quantitative decay 
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properties of solutions of a linear wave equation and led to a small data global existence in 
dimension n > 3 and an almost global existence in R 3+1 . Furthermore, in IR 3+1 , a sufficient 
condition on a quadratic nonlinearity F, which guarantees the small data global existence, is 
the celebrated null condition introduced by S. Klaincrman [14]. Under this condition, D. 
Christodoulou [3] and S. Klainerman [16] independently proved the small data global existence 
result. 

The approach of [3] used the conformal method, which relies on the embedding of Minkowski 
space to the Einstein cylinder Rx S 3 . In [16], S. Klainerman used the vector field method based 
on the weighted energy inequalities generated by the vector fields 

T = {Qij = x t dj - xjdi, Li = tdi + Xid t ,d a ,K = (t 2 + r 2 )d t + 2trd r , S = td t + rd r ) ■ (4) 

The vector field method found many applications, in particular to systems of nonlinear wave 
equations on R 3+1 with multiple speeds [18], [27], [28] and obstacle problems [22], [23]. Due 
to the lack of Lorentz invariance, these works avoided the use of the hyperbolic rotations Li, but 
still used the scaling vector field S. 

Another application of the vector field method is to the quasilinear wave equation of the form 

n gW </> = F(4>,d<f>) (5) 

with g(0) = to, the Minkowski metric. The quasilinear part of the equation g a/3 (cf))d a p<f> never 
satisfies the null condition. Nevertheless, several examples of (5) still admit global solutions. In 
[19], [20], H. Lindblad obtained the small data global existence result of the equation 

g a ^<f,)d a p4> = 

on R 3+1 . A particular case 

d tt 4>-{i + (t>) 2 ^ = {) 

was investigated previously by S. Alinhac [1]. 

Separately, the motivation for studying (5) arises from the problem of global nonlinear stabil- 
ity of Minkowski space originally proved by Christodoulou-Klainerman by recasting the problem 
as a system of Bianchi equations for the curvature tensor [4]. Later, Lindblad- Rodnianski [?] 
obtained a different proof of stability of Minkowski in wave coordinates, in which the problem 
was formulated as a system of quasilinear wave equations for the components of the metric. 

We should also mention that a linear problem 

U g(t,x)4> = 

has been studied in [2]. There it was shown that <j> has the decay properties similar to those 
of a solution of a linear wave equation on Minkowski space provided that g(t, x) approaches the 
Minkowski metric to suitably as t — > oo . 

A common feature of these problems is that the background metric g converges to the 
Minkowski metric either by the assumptions in [2] or, for the equation (5), by the assumption 
g(0) = m and the expected convergence </>(t, x) — > as t — > oo. The need for such convergence, 
or at least convergence to some stationary metric g(x) is dictated by the vector field method. All 
applications of the vector field method require commutations with generators of the symmetries 
of Minkowski space, at the very least with the scaling vector field S = td t + rd r . For the problem 

the error term coming from the commutation with S (or Li) would be of the form tdtg a ^d a p<j) 
which leads to the requirement that tdtg is at least bounded and thus g decays to a stationary 
metric. 
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Therefore, the vector field method has not been so far successful in application to the equation 



u g(t,x)^ = F{<t>,d4>), ^ 
0(0, x) = e(f> o {x),d t (j)(0,x) = e<t>i(x), 

in which g(t,x) does not converge to a stationary metric g(x). These problems describe propa- 
gation of nonlinear waves in a time-dependent inhomogeneous medium. 

In this paper, we develop a new approach to treat the small data global existence problem for 
the equation (6). This approach relies on a new method for proving decay for linear problems, 
developed by M. Dafermos and I. Rodnianski in [5] . We now describe the assumptions and the 
main results. 



We assume that the metric g coincides with the Minkowski metric outside a cylinder 

g a p=m aP + h aP , g afi = m a ? + h a(i , 

where h a p are smooth functions supported in {(t, x)||x| < \R}, R is a fixed constant and m is 
the Minkowski metric. We make a convention that the Greek indices run from to 3 while the 
Latin indices run from 1 to 3. 

In coordinates (t, x), we denote do = dt, 9, = d Xi , 8 = (81,81,82,83). We also use the 
standard polar coordinates in Minkowski space (t, r, w). Let denote the induced covariant 
derivative and ^ the induced Laplacian on the spheres of constant r. And Q is a shorthand for 
the angular momentum = Xidj — Xjdi. We also define the null coordinates u = v = ^jp. 
And denote T as the vector field 8 t in (t, r, w) coordinates. 

In our argument, we estimate the decay of the solution with respect to the foliation E T , 
defined as follows: 

S T := {u = u T , v > v T }, 
S T :={t = T,r<R}L)S T , 

where u T — v T — Thus the corresponding energy flux in Minkowski space is 

E[4>]{t) := / \8 t 4>\ 2 + \8 r ^\ 2 + \f<j>\ 2 dx + [ (\8 V 0\ 2 + \f0\ 2 ) r 2 dvdu:. 

Jr<R JS T 

Additional to the assumption that g coincides with the Minkowski metric outside the cylinder 
{(t, x)\\x\ < \R}, we make two other assumptions: 

Al : There exists a positive constant A such that in (t, x) coordinates 

goo < - A, A|x| 2 < gijXiXj < A _1 |x| 2 
for any x = (xi, X2, X3) € K 3 . 

A2 : There exists a positive number a < 1 and a constant Co such that two integrated local 
energy inequalities hold for any smooth function <fi(t, x) 



'Ti Jr<^R 

where we denote 



/ 2 / |90| 2 + —dxdt < CoEWfa) + CoD[n g <f>]%, (7) 

Jti Jr<R r 

l^l^xd^Co^^^rO+^^n^^+SM^O+^p^]^), (8) 



D[F] T T \:= f 2 f \F\ 2 (l + r) a+1 dxdT. 
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We define two sets 

A:={(k,j)\k + j<8,k<5}, 
B ■= {(k,j)\(k,j + 2) e A}, 

where k, j are always nonnegative integers. And then denote 

(k,j)eA 

which is determined by the initial data <j>o(x), (f>i(x) and the equation (6). 
Our main results are: 

Theorem 1. Suppose the nonlinearity F satisfies the null condition (2) and g satisfies Al and 
A2. Assume that the initial data (f>o(x), 4>i(x) are smooth and supported in {\x\ < R}. Then 
there exists eo > 0, depending on R, a, Eq, X, Co, h, such that for all positive e < eo, the 
equation (6) admits a unique global smooth solution. Moreover, for the solution <f>, we have 

(1) Energy decay 

E[<f,](T) <CE e 2 (l + T)- 2+a 
for some constant C depending on R, a, X, Co and h. 

(2) Pointwise decay: for any a < S < 1 

\<t>\ < C s ^e(l + r)-\l + \t-r + R\)~i+i s , 
J2 < C^Eo~e(l + r)"5(l + \t-r + R\)- 1+ ? a , 

l/3|<2 

where Cs also depends on 5. 

Remark 1. The argument here is also applicable to the corresponding problems in higher di- 
mension without null condition. 

Remark 2. The same conclusion holds if the assumption A2 is replaced with that g satisfies (7) 
and the deformation tensor irj^ = dtg a p is small, independent of the initial data. We remark 
that this is consistent with our attempt to investigate nonlinear wave equations on backgrounds 
far from Minkowski space. 

Remark 3. It is not necessary to require that the initial data have compact support. The general 
assumption on the initial data can be that the following quantity 

V // r 2 - a \dfl k T3<j)(0,x)\ 2 dx 

(k, 3 )eA JJ ^ 

is sufficiently small, where a comes from Ai or the smallness assumption (19). 

The condition Al is easy to satisfy. It is equivalent to say that the background is uniformly 
hyperbolic and the vector field dt is uniformly timelike. Therefore we reduce the small data 
global existence to the two integrated local energy inequalities. Below, we describe particular 
conditions for which Ai can be explicitly verified. We will show that under the assumption that 
g is C 1 (M 3+1 ) close to the Minkowski metric, A\ and A2 hold and hence Theorem 1 follows. 
More precisely, denoting 

H = m&x{\\h al3 \\ c i ,\\h al3 \\ c i} , 

then 

Theorem 2. Suppose h is supported in the cylinder {(t, x)\\x\ < -j}. Then there exists a positive 
constant eo, depending only on R, such that if H < eo, then g satisfies conditions Al and A2. 
Hence Theorem 1 holds for some a < 1. 



4 



Remark 4. From Theorem 2, we retrieve the classical result proved by Klainerman [16] and 
Christodoulou [3] in Minkowski space where H = 0. 

The main difficulty of considering nonlinear wave equations in inhomogeneous media is the 
lack of symmetries compared to Minkowski space. In order to make use of those symmetries 
in Minkowski space, previous works have relied on the fact that the background metric decays 
to its stationary state, which is not satisfied in this context. In fact, we even allow the back- 
ground metric g to stay far from the flat one provided that we have two integrated local energy 
inequalities. 

In our approach, we avoid the use of vector fields containing positive weights in t, e.g, 
S = tdt + rd r , Li = Xid t + tdi. Traditionally, vector fields from the set T (4) are used both as 
multipliers and commutators. In this paper, we only commute with d t , flij- The role of multiplier 
vector fields is played by d t , /J^, where / is some appropriate function, in the derivation of the 
integrated local energy decay, and by the family of vector fields r p (d t + d r ) localized to a far 
away region r > R. The key in our argument is a new approach, developed by M. Dafermos 
and I. Rodnianski in [5], to the problem of decay, in particular, of the energy flux E[4>](t) for 
solutions of linear wave equations. This new approach is a combination of an integrated local 
energy inequality and a p-weighted energy inequality in a neighborhood of the null infinity. We 
will discuss them in details in Section 2 and Section 3. 

The plan of this paper is as follows: we first establish an integrated energy inequality in the 
whole space time in Section 2 by using the vector field method; then in Section 3, we revisit the 
p-weighted energy inequality developed in [5] and prove the decay of the energy flux £?[</>] (t). 
In Section 4, we use elliptic estimates to get the pointwise decay of the solution; and then in the 
last two sections, we close our boostrap argument and conclude our main theorems. 

2 Integrated Local Energy Inequality 

In this section, we use the multiplier method to prove an integrated energy estimate in the whole 
space time under the conditions .41 and .42 or the smallncss assumption on H. This kind of 
estimates were first proven in [25]. Here we would follow the way in [6], also see [30]. 
We recall the energy-momentum tensor 



We denote n as the unit normal vector field to hypersurfaces and da the induced measure, 
dvol the volume form of (R 3+1 ,g). Denote the null infinity from n to T2 as 




Given a vector field X, we define the currents 



where tt* = \ 




X T T \ := {(u,v,u})\u Tl < u < u T2 ,v = oo} 



and the corresponding energy flux 




v = 



-oo 



which can be interpreted as a limit when v — > oo. 



•5 



Define the modified energy 

E[<t>]{r) = E[<t>]{r) + 
Take a vector field defined as follows 

X = fd r = f^-di, 

r 



where / is a function of r. Consider the region bounded by the hypersurfaces E Tl and S T2 . 
Using Stake's formula, we have 

f jx^da- ( jx^da- f J?[<t>]n»dv 
= f 2 f D»J x [4>]dvo\ = f 2 f FX{4>)+K x [cj)]dvo\ (9) 

where 

K x [<f>] - V"[<j>]ir x v - <)yJ±-uVo ■ - (if + r-\f)d^d^ 

+ h-dig^ ■ d^d v 4> ~ -J-d i9llv ■ gr&w^, 
2 r 4 r 

in which we denote /' as d r f. 

Choose another radial symmetric function \ of r. We have the following equality 

-xd^d^ + ±n gX ■ 4? = -i x (n^ 2 - 2<j>u g <j>) + ±n gX ■ 4> 2 
= \ ( a 9 x ■ 4> 2 - x^ 9 4> 2 ) + X^g4> 

Add the above equality to both sides of (9) and modify the current as 

J?[4>] = J*[4>\ - \d»x ■ + \xd^ 2 - (10) 

Then we obtain 

f J x [4>]n»da- f jx^da- f T J*[cj>\n»dcj (11) 
= jT £ FX(<f>) + Fct> X +(x- r-\f + \mdAf -\u gX - <f 
+ (r-\f + \f - X)(dt& 2 + (X- \f')\n\ 2 + error, 



where 

error = dj(f^)h^d^ ■ drf - {\f + r" 1 / - xW ^ + |5^| 2 - \0^\ 2 ) 

- \i-di<T ■ M^<P - -J-d i9liV ■ gTff 1 ^. (12) 
2 r 4 r 

Here recall that h af} = g a ' 3 - . 

The idea is that we choose functions / and x such that the coefficients on the right hand 
side of (11) are positive. And then control the left hand side by the energy E[4>]. Thus we end 
up with an integrated energy inequality in the whole space time. To proceed, let's first prove 
several lemmas in order to estimate the left hand side of (11). 

We first show that the spherical average of the solution near the null infinity can be bounded 
by the energy. 
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= 0. 

V — oo 



Lemma 1. Let <j>(t,x) G C°°(1R 3+1 ). Then we have 

rf \4>(t,r 7 uj)\ 2 duj < E[4>](t), V(t, r, w) e S T U I r . 

Proof. Suppose the point can be represented as (u, v, uj) = (t, r, uj) in our coordinate systems. It 
suffices to consider the case when E(t) is finite. We claim that <j) vanishes at the null infinity 
Iq t . In fact since I[4>]q is finite, we conclude that J x u T <\>\da is finite. Recall that r = v — u. We 
infer that 

/ ( / \d u 4>\du\ duj< [ <p 2 u do- ■ [ r- 2 du 

Hence d u <j) vanishes on X% T . Notice that the initial data are supported in r < R. Then the finite 
speed of propagation for wave equations [29] implies that <p vanishes on So- We thus conclude 
that <j> vanishes on Xq t . Therefore 

J \<p\ 2 duj = J ^jf dvtpdvj duj<J 4>l<to-£ r- 2 dv<^-E[4>]{T), (13) 
where on S T , da = r 2 dvduj. Hence the lemma holds. □ 

The following analogue of Hardy's inequality will be used later on. We borrow the method 
from [6] . 



Lemma 2. If (p is a smooth, then 

2 n / j \ 2 



(14) 



In particular 

[ ct) 2 dx < 6(1 + R) 2 E[(t)}(T). (15) 

Jr<R 



lr<R 

Proof. Take a function rj as follows 



rj(r) = r - 21n(l + r) + — 

1 



+ r 
Then 

r?/ W = (TT7)2' ^°) = ' l»?( r )l<»-. 

Denote dcr as dx when r < R and r 2 dvduj when r > R. Integration by parts and using lemma 
1 imply that 



2 /■ /-Ji 



dcr = I I (f) 2 dr]du) + (j) 2 drjduj 



oo 



= / <j> 2 nduj — 2 / 4> r drdw — 2 / ryi/i • 4> v dvduj 

<%](r) + J/ n 2 r-^ 2 da + 2 f $.dx + l[ <gdo 
<3%](r) + ^ n\~^ 2 da, (16) 



where r — v — u and on S T , u is constant. 

Notice that n is nonnegative and ln(l + r) > j^. We conclude that 

n _ r / hi(l + r) 1_ \ < r 



r 1 + r V r 1 + r / 1 + r 
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Hence 



da 



The inequality (14) then follows from (16) by absorbing the second term. Inequality (15) 
follows from (14) if we restrict the integral in the region r < R. □ 

In the region r > R, we analyze weighted solution ip = r<j> instead of <j> itself. In the energy 
level, they are equivalent in the sense of the following corollary. 

Corollary 1. In the outer region r > R, we have 

[ \d v (r<t>)\ 2 + \y{r^)\ 2 dvdu:- ( ((d v <j)) 2 + \f^\ 2 ) r 2 dvdcj 
Proof. In fact, since 

/ (r4>)ldvdu= [ [r 2 <j) 2 v + {r(b 2 ) v ]dvdLO= ( <gda + / rtf 
the corollary follows from Lemma 1. □ 



< 2E[</>](t). (17) 



Now we are able to estimate the boundary term in (11). 
Proposition 1. Suppose f and \ satisfy 

\f\<cu \x\<£- r , \x'\<j^ 



for some constant C\ . Then 



^da 



<6d{l + H + 2H 2 fE[<j)){T). 



Proof. For r < R, notice that on the surface t = r, we have 
jX[<j>]n»da = -jx^WV^Gdx 

Since \g a0 \ < 1 + H and | 5 a ^| < H, a ^ (3, we have 



-Gdx. 



< a/(1 + H + 2H 2 ) 4 = (1 + H + 2H 2 ) 2 



and 



<f>d r <f>\ = \fg^d^d r ^\ < l -±^-C 1 \dcj > \ 2 . 



Under the assumptions on \, we can estimate 



and 



\\x -dV\ = Ix^M <gi 1+ 2 4g ( if r)a 



Ci 



(1 + r) 2 
1 + H, nil2 



(18) 



For r > R, the metric is flat according to our assumption. The unit normal vector field to S T is 

n=-L dv = -L(d t + d r ). 
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Thus we can calculate 



jX[<t>]n»da = Q/(M 2 - |y0| 2 ) - \ X ' • 4? + X ■ d v cj> ■ ^ r 2 dvdw. 



On the other hand, the condition on \ shows that 



Therefore, according to (14), we can show 

/ J* [(fin* 1 da <d(l + 3H) [ \dcf>\ 2 V^Gdx + ^-I / (|9,,0| 2 + |^| 2 ) r 2 d«(iw 

< ^1(1 + H + 2H 2 ) 4 E[4>}(t) + ^1(1 + H + 2H 2 ) 4 ■ 6E[4>](t) 

< 6d(l + H + 2H 2 ) a E[4>](t). 

Hence we conclude the proposition. □ 

Remark 5. If E[<P](t) is finite, all the above statements hold if we replace E[^)\(t) with E[4>](t) . 
The reason is as follows: under the assumption that E[<p]{r) is finite, Lemma 1 holds for E[4>](t) . 
Hence Lemma 2, Corollary 1 and Proposition 1 also hold if we replace E[4>](t) with E[4>](t). 

We hope that, by using the multiplier method, we can derive the H 1 estimates of the solution 
on E T . That is we want to show that the energy current ,/J[0]n A ' has a positive sign. The 
proposition below guarantees this positive sign provided that (fir'- 7 ) is uniformly elliptic. In 
addition, the elliptic method we will use later on also requires the uniform ellipticity of (ff*- 7 '). 
This is obvious when (!R 3+1 , g) is a small perturbation of Minkowski. We claim that it still holds 
under the condition Al. 

Proposition 2. If g satisfies Al, then there is a constant Ai such that 

Ai < -9 00 < Ar\ Ax/axg < (g ij ) < X^hxs, 

where (g^ u ) — (g^)" 1 and u, v run from to 3 and i,j run from 1 to 3. 

From the geometric point of view, the hypersurface t = const is spacelike. Therefore its 
normal is timelike. Since we require d t to be timelike, the positivity of the current jJ[0]n M 
follows from the fact that T(X,Y) is positive for any two timelike vector fields X, Y. However, 
we are concerned about the uniform lower bound. We prefer the following algebraic proof of 
Proposition 2. 



Proof. We write the matrix g as 



~a b 
b T D 



where a = -g 00 , b = (301,502,503) and D = (gij). Then G = det(g^) = -(a + bD V )detD 
and 



or) 



, \ / 1 bD~ 

-a b 



a+bD- 1 b T a+bD- 1 b T 
b T D ~\ D- 1 b T n-1 _ D- 1 b T bD- 1 

\ a+bD- 1 b T ^ a+bD- 1 b T 



Since g satisfies condition Al, we have —a = —.goo > A and D is positive definite. Thus 

A< a + bD~H T < l + iJ + A^H/jH 2 <1 + H + 'i\- 1 H 2 . 
Denote X = bD~ x . For any vector Y, we have two cases: 
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If Y is perpendicular to X, then 



D- 1 h T hD- ls 



IfY = X, then 



, i D~ 1 b T bD~ 1 \ T aXB- 1 X T XD~ 1 X T XDX T - \X\ 4 



a + bD-W J a + bD-W a + bD' 1 ^ 

Notice that 

XD~ 1 X T XDX T > \X\ 4 , 

which follows from Cauchy's inequality if we assume, without loss of generality, that D is 
diagonal. Therefore we have 

AW <r fir' -£™£V 



l + H + SX^H 2 ~ V a + bD~ l b T 



Now g 00 = — a+h p-i b T and (g lj ) = D 1 — ^p^-rrgr- , the Proposition then follows if we choose 



\ - x 

A1 l+H+ZH'^X- 1 ■ 



□ 



Having proven this proposition, we conclude that the energy flux through E T in (R 3+1 , g) is 
equivalent to that in Minkowski space. 

Corollary 2. If g satisfies condition Al, then there is a constant A 2 such that 

A 2 £[0](t) < f jJ[0Kda < X^EMr). 
Proof. In (t, x) coordinates, on {|a;| < R} n S r , we have 

= \ {d l <j>di4> - d^dtcf)) V^Gdx 

Proposition 2 then implies that 

3 3 
Ai ]T |9^| 2 < jidiWrf ~ 9 m %4> < A- 1 \M\ 2 - 

(8=0 0=0 

On one hand, using the notations in Proposition 2, we have the lower bound of — G 

-G=(a + bD- 1 b T )detD > A 4 . 
On the other hand, we have already shown in (18) that 

|G| < (1 + H + 2H 2 ) 4 . 

A A 2 

Let A2 — 2 (i+3H) 2 - ^ e °ktain 

p 3 /* /* 3 

A 2 / Y^\d^\ 2 dx< Jl\(j>]n^du<\^ Y,\ d ^\ 2dx - 

Jr<Rp =Q Jr<R Jr<Rp =0 

Notice that the metric is flat on S T . Thus we have 

jJ[^Kda = l - (\d v ^\ 2 + \f^\ 2 ) r 2 dvdu. 
Since A 2 < i, the Corollary then follows. □ 
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The proof of Theorem 1 is quite similar to that of Theorem 2. We thus prove them together. 
For Theorem 2, we prove it by taking 



eo = sup 



<i 700(1 + \R)W 



2 

Thus for H < e , there exists a positive constant a < 1 such that 

a 

H ~ 700(1 + \R) a+1 ' ^ 
From now on, we instead use condition (19) on H. 

Having the above basic preparations, we are now able to establish the key estimates in this 
paper. 

Proposition 3. If g satisfies the conditions Al and A2 or the smallness assumption (19), then 
there is a constant C depending on H, i?, a, X, Co such that 



(1) Integrated energy bound 

fT2 r ^2 _,_ \p, ^|2 _l iv", ,|2 
In JT, T 

In particular, for the smallness assumption (19), C can be 

(2) Energy bound 

%](r 2 ) + /[C? < C%](ri) + C£>re. (21) 

(3) Improved integrated energy bound in angular directions 

f 2 j ^-dxdr < CE[<t>]( n ) + CD[F}1\. (22) 

Jti JS t r 

Here recall that D [F] £J = \F\ 2 (l + r) a+1 dxdr. 

We mention here that variants and generalizations of estimate (20) can also be found in [?], 
[23]. To be consistent, we give the proof of the above proposition by following the method in [6]. 

Proof. We choose the function /, \ suitably to make the coefficients on the right hand side of 
(11) positive. Take a exactly as the one in condition A2 or (19) and f3 = — . Set 

f = @ ~ (2 + r )a ' X = r /• 

Direct calculations show that 

r-'f+h'-X ' 



x' = -P 



Y (l + r) ^ 1 ' 

(l + r) a+1 -(a + l)r-l 



r 2(l + r )a+l 



x// = 2 ^ r -3_ 2/3 4 2(« + l) 



r 3 (l + r)« r 2 (l+r) a+1 r(l + r) 
2(a + l) x 4 xj 



a+2 - 



vX r(l + r)«+ 2 r 2 X \ r 



Notice that 

< " " < (1+0 



a / (i + r-r-i ^ , 



1 + r 
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We obtain 



(3 l + r(l + r) a -l (3 
X=-r—— M , < 



1 + r r (1 + r) a ~ 1 + r ' 

X - r f+nf 



a+l ■ 



2 J (l + r)^ 1 (1 + r) 
v V_ r -i* gO + O+I _ /3((l + r)"-l) 1 1 



2 r(l + r) Q+1 r(l + r) a (1 + r) 1+tt ~ (1 + r) Q+1 ' 

To apply Proposition 1, we need to estimate X ' . Since 

Q ^ (l + r )°+i-( a + l)r-l = (l + r)»-l + (l + rT-or-l ^ ^ + r)a _^ 



we conclude that 



\X'\ = -X'< n -(l + r)°- J - 



(l +r )a+l ^ ' (1 + r) 2 ' 

Notice that |/| < /?. We have shown that /, x satisfy conditions in Proposition 1. 

We have two cases according to the conditions g satisfies. 
Casel : If h satisfies the smallness condition (19), then 



-G 

= -9 ij d ijX - (d a g m + \g m d a g M ■ g^ d lX 

> -A X - H\d ijX \ - {AH + 32ff(l + Hf) ^| x '|. 



As having computed above, we find that 



2 , 2(a + l) 
-A X = -x" - -X - 



r(l+r)«+ 2 ' 

^ r (/ + r)«+ 2 + ' X ''r " 7(1+7F+^' 

Recall that H < 700 ( 1+ i^^ + i , a < 1 and h is only supported in r < \R. We conclude 
that 

1 

~ D9X - r(l +r)«+ 2 ' 

It remains to show that the error term (12) can be absorbed provided that H is small. In 
fact since 

' ' ~~ (1 + r) a+1 1 /l_ (l + r)« r S (1 + r)» S ' 

we have 

Mf*) = \f'^ + ffe- E § L )\<2- 

Using Cauchy-Schwartz inequality, we obtain 

E ^ £(i^i 2 + \\ d ^ + \ £ i^i 2 + i^i 2 ^ f 150,2 

i i,j 

and 



19^0+ \d t <j>\ 2 - |^«/>| 2 | < l^^^^l < 4|50| 2 . 
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Therefore we can estimate 



\error\ < 6H ■ ^-\d(j>\ 2 + 4H\dcj)\ 2 + 4/3H\d4>\ 2 + 2/311(1 + 4H) 2 \d(j)\ 2 
< l8/3H\d^\ 2 < ni+ \ R)a+ ^ d ^ 2 + + m\ 2 ) 



< 



1 



(( dt <t>y + {d r <t>y + \nn 



18(1 + r) c 

where we recall that that f3 > 2, H < 700 ( 1+ 1 j^c+i an d the error is supported in r < \R. 
On the null infinity , the unit normal is 

n=±d v = ±(d t -d r ). 

Thus we can write 

J*[4>Kdo- = \ (f(-(d u 4>) 2 + m 2 ) -x'-^ + Zx- d u <j> ■ tf>) r 2 dudu. 
Notice that at null infinity, x' — and / = (3. The above inequality implies that 



F ° r H ^ 700(l4^)° + 1 ' ^ haVe 



-G > (1 - H) 4 - 6(1 + HfH 2 - 8(1 + H)H 3 -9H 4 >(l-H- 2H 2 ) 4 . (23) 
Using Proposition 1 and the inequality (11), we obtain 

J T1 ^ (l + r)«+! r(l+r)«+ 2 V ' 

< 7(3 (%](n) + %](r 2 ) + ^[C?) + 2/3 jT jf |F«M + iF^-ldxdT. 

Case2 : If g satisfies A2, although we can not guarantee the positivity of — D g x as in case 1, we 
can control this term and the error term by using the condition (7). The error term can be 
estimated in the same way as having shown in case 1. We thus obtain a similar inequality 
to (24) but with a constant multiple of _D[F]pJ on the right hand side . Although this 
constant depends on a, R, A, h, Co, the explicit constants are not essential for this case. 

Now, we need to relate the energy flux on E Tl and that on S T2 . For this purpose, take X to be 
the vector field T. We have 

\K T [<j>}\ = \ld t g^T^\ < y (4 + 2(1 + 4H) 2 )|^| 2 < 3H(1 + H + 2H 2 f\8<P\ 2 . 

By (9), we obtain 

/ jJ[</>Kdcr+ f jT[ ( j ) ]n IJ -do-= [ jj[0]n"da - I 2 f F<j> t + K T [frivol 

JZ T2 Jill "' S -l JT, t 

<[ J^[4>}n^da+ [ 2 f \F(/) t \dvol + 3H(l + H + 2H 2 ) 5 [ * [ \d4>\ 2 dxd T . (25) 



Notice that 



J^[(p]n^da = i (d^dicj) - 0*05^) v^Gdx. (26) 



Again we distinguish two cases. 
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Casel :When H satisfies (19), we can show 

|0W - - (d^f\ = fruloiho - g 00 (d^) 2 - |^| 2 | 

< H\8cj)\ 2 + 2H(d l( j)f < 3H\d(t)\ 2 . 

Using the facts that the metric is fiat outside the cylinder {(£, x)| < jR} and H < 
inequalities (18), (23) and (26) imply that 

^E[<f>](r)<2 [ J>]n^<^](r). 



Recall that E[4>]{t) = E[4>](t) + I\<f\l. By (25), we have 

1 7 1 3 f T2 f 20 r T2 f 

e[4>](t2) + jg/tc < g^M(ri) + - y iF^idcdr + y F y } <x w?**r. 

Restricting the integral of the left hand side of (24) to the region r < \R,we can estimate 
the above inequality as follows 

7 ~ 13 f T2 f 20 a / - - \ 

< gEW(n) + y J ^ + -—TP (EMM + mM) 

which implies that 

%](r 2 ) + < 2%](n) + lf 2 J 



Together with (24), we can show that 



i ri is T (l + r)«+! r(l+r)«+ 2 

< 21/3%](n) + |s jf" jf 7|F&| + |F0 r | + IF^-Idxdr. 



/ti 

Using Cauchy-Schwartz inequality to |-F0 r |, |-FV>t|, ^i+f > f° r example 



2 2 y 2 y 2 (1 +r) a + 1 

we can show that 



6? + 02 + 



(l + r )a+l r (l +r )«+2 



dxdr 



(27) 



< 42/3%](r 1 ) + (^f/?) 2 + 2(3/3)^ £>[fl£ (28) 

<42/3^[0](r 1 ) + (18/3) 2 J D[0]^. 

Hence (20) holds with constant C — 18 2 /3 2 for this case. This explicit constant will be 
used later on. The energy inequality (21) follows from (27), (28). 
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Case2 : When g satisfies conditions Al and A2, condition (7) together with Corollary 2 and 
estimates (25) lead to a similar inequality of (27) but with a constant multiple of D[F]^? 
on the right hand side. This constant also depends on Co (in A2), H, A2(in Corollary 2). 
However, once we have shown (24) and (27), the estimates (20), (21) follow similarly as 
in case 1. 

We hence have proven (20), (21) in any case. Finally, the estimate (22) follows from the 
fact that 

!,/_ -io (3(l + r)+r R 1 

X 2 J T P r(l + r)«+! " (l + i?)«+! r' T ~ 

□ 

It will be much more convenient if we use E[4>] (r) instead of E[4>] (r) in the following argument. 
We use E[4>](t) in order to make the previous argument rigorous. However, as we have pointed 
out in Remark 1, all the statements above hold if we replace E{<P](t) with E[4>](t) provided that 
E[4>](t) is finite. The idea is that under the bootstrap assumptions, we first show that E[4>](t) 
is finite. And then we conclude that Proposition 3 holds if we replace E[4>](t) with E[4>](t). 

Corollary 3. Proposition 3 holds if we replace E[4>](t) with E[4>](t). 

Proof. If D[F]q is finite, then E[c/)](t) is finite according to (21). Thus by Remark 1, we infer 
that Proposition 3 is also true if we replace E[4>](t) with E[4>](t). 

If D[F]l is infinite for some t, then Propositiuon 3 holds automatically if we replace E[c/>](t) 
with E[4>](t). ' □ 

Restricting the left hand side of (20) to the region r < R, we obtain an integrated local 
energy inequality. 

Corollary 4. If g satisfies (19), then 

A2 



'1T 2 
J-Tl ' 



(29) 



/ 2 / $ + \d r (f>\ 2 + l^l 2 + — dxdt < CE[$\{ Tl ) + CD[F] 

J Tl Jr<R r 

where C = 18 2 /3 2 (1 + R) a+2 , /3 = f . In particular 

f 2 f $ + \d r cf>\ 2 + \f^\ 2 dxdt < CzElflin) + C 2 D[F}H (30) 

Jti Jr<R 

with C* 2 = 18 2 (3 2 (1 + R) a+1 . 

If g satisfies (19), we can show that g satisfies the condition Al. The corollary implies that 
g also satisfies the condition (7). 



3 Weighted Energy Inequality 

In this section, we revisit the p- weighted energy inequality in a neighborhood of the null infinity 
originally developed by M. Dafermos and I. Rodnianski in [5]. Starting from the integrated 
local energy inequality (29), the p- weighted energy inequalities allow us to obtain the decay of 
the energy flux E[4>](t). Moreover, the p-weightcd energy inequailites play an important role in 
estimating the quadratic nonlinearity with null condition. 

To avoid too many constants, we use the notation A < B to indicate that there is a constant 
C, depending on R, a, A, h and Co, such that A < CB. 

Proposition 4. For any < p < 2, we have the p-weighted energy inequality 
f r p ipldvdu + [ 2 f r 1 " 1 (iPt + lf iP\ 2 )dvdcudT 

JS T2 J Tl JS T 

<E[cj ) }{T 1 )+ f rP^dvdLo + [ 2 [ rP +1 \F\ 2 dvol + D[F} T T l, $ = r<j>, (31) 
where the implicit constant also depends on p. 
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Proof. When r > R, the metric is flat. Rewrite the equation (1) in null coordinates 

-d u d v ip + fcip = rF, V := r<f>. (32) 

Multiplying the equation by r p d v ip and then integration by parts in the region bounded by the 
two null hypersurfaces S Tl ,S T2 and the hypersurface {(t, x)\\x\ = R} } we obtain 

/ r p (d v ip) 2 dvduj + f 2 f 2r p+1 F ■ d v ipdvdrduj 

JS T2 Jn JS t 

+ J 2 J r^ 1 (p(d v ^) 2 + (2 -p)\fi>\ 2 ) dvdrduj + J ^ r p \ftP\ 2 duduj 

= f r p (d v ^) 2 dvduj + I ' r p (|y^| 2 - (d v ^) 2 ) du:dT\ r = R . 
JS T1 Jn 

We claim that we can estimate the boundary term on {r = R} as follows 

r (|yV| 2 - (dvipf) dujdr 

J Ti 



(33) 



(iy^| 2 - (d v ^) 2 ) dud T \ r=R 



R p 



< 



£>re+£7[0](n). 



In fact, it sufflces to consider the case for p = 0. Thus set p = in the above p- weighted energy 
inequality (33). Notice that 



[ T2 f 2r\ 

Jn JSt 



F\ ■ \d v tp\dvdTduj < 



\F\ 2 (l + r) 



ct+l 



'n JS T 
<E[4>](n)+D[F] 



, r 2 tf,+<p 2 
J r 2 {l + r) a + 1 



dvo\ 



and 



/ rP\f^\ 2 dudu= ( \f(f>\ 2 r 2 duduj </[0]; 2 . 



Inequalities (17), (21) and (22) then yield the desired estimates. Here we recall that, by Corollary 
3, Proposition 3 holds if we replace E[4>]{t) with E[<fr](r). 

For general < p < 2, we can estimate the inhomogeneous term as follows 

r p+1 \Fd v iP\ < V- 1 ^) 2 + -r p+3 \F 2 \. 
2 p 



Then the p-weighted energy inequality (31) follows by absorbing the first term. 



□ 



Intuitively, the inequality (31) indicates that the solution 0, at least the good derivative 
of the solution d v (j), decays fast in r since p e [0, 2) and ip — rcj). However, for the nonlinear 
problem, merely decay in r is not enough to obtain the long time existence of the solution. Wc 
thus have to transfer the decay in r to the decay in time t. The first step to realize this is to 
find some quantities which decay in t. However, in this context, decay in time t is characterized 
by the parameter r of the foliation E T . We show that, under appropriate assumptions on the 
nonlinearity F, the energy flux through the hypersurface S r decays in r by combining the 
integrated local energy inequality and the p-weighted energy inequality. 

Proposition 5. If there is a constant C\ such that F satisfies the following conditions: 

(a) /^/ ST |F|V-«»dvol<C 1> 

(b) J^J^lF^dvolKC^l + n)-^, 

(c) D[F] T T \ < Ci(l + n)~ 2+a + E[4>](ti) 
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for all T\ < t 2 , then we have the energy flux decay 

E[d>}(r)< (e 2 Eb + C 1 )(l + r)- 2+Q . 

Proof. We first take p = 2 — a in the p- weighted energy inequality (31). Since the initial data 
are supported in {\x\ < R}, the finite speed of propagation for wave equations shows that the 
solution <f) vanishes on S - Let r 2 = r, n =0. Under our assumptions on F, we obtain 



/ r 2 ~<*tfc 



dcudv < e 2 E +Ci, (34) 



/ 2 / r^ip^dcodvdT < [ * [ r^ip^dujdvdT < e 2 E + C x . 
Jt x Js t Jo Js t 



(35) 



The second inequality holds for all t\ < t 2 . We claim that there exists a dyadic sequence {t„}^L 3 
such that 

r^^dvdio < (1 + r n )- 1 (e 2 E + Ci) , (36) 

where r„ satisfies the inequality 7~ 2 t„ < t„_i < 7 2 t„ for some large constant 7 depending on 
R, A, a, /i and Co- In fact it suffices to show that there exists r„ <E [7™,7™ +1 ] such that (36) 
holds. Otherwise 



/ 7 fc Js. 

which contradicts to (35) if 7 is large enough. 
Interpolate between (34) and (36). We get 



I 7 I r 1 - a i; 2 v dvdu J dT>\n 1 (e 2 Eo + C 1 ), 
J~i k Js^ 



f ri> 2 v dvdu<(l + T n )- 1+a (e 2 E + Ci). 



(37) 



Now, take p = 1 in the p-weighted energy inequality (31). Using the estimates (37) and 
conditions (6), (c), we obtain 



/ r{d v ip) 2 dujdv + I [ (d v ip) 2 + \fip\ 2 dujdvdt 

J S't J T n -i J St 

< f r(a^) 2 d W d« + J BM(r„_ 1 ) + (l + r„_ 1 )- 1 + a C 1 

< (1 + r n _!)- 1+a (e 2 E + Ci) + S[0](t„_i) 

< (1 + r)- 1+Q (e 2 E + d) + £?[^](r n _i) 

for all r e [t„_i,t„]. 

Our goal is to retrieve the full energy flux through the hypersurface E T . Since 

r<p 2 duj 



(38) 



{d v ^) 2 + \fip\ 2 dudv = [ (<p 2 v + \f<f)\ 2 )r 2 dLodv- 

> 2 f Jl[<t>]n»da-R [ 4> 2 du 

and 

rR rR rR. rR. rR. 

R 2 <j) 2 =2 r(/) 2 dr + 2 r 2 (f> ■ d r (f>dr < 2 / r(/) 2 dr+ r 2 4> 2 dr + r 2 (d r (f>) 2 dr, (39) 
Jo Jo Jo Jo Jo 

integrate over [ti,T2] x S 2 . We obtain 

R 2 f [ (j) 2 (T,R,u})du}dT < f f 2 — +4> 2 + (d r (f>) 2 dxdT 

J Tl J\u\ = 1 Jtx Jr<R r 

<£[0]( ri ) + (l + ri )- 2 +«C' 1 
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by (29) and condition (c). Adding (29) and (38), we can show that 

f " E[</>](T)dT = f " f 4> 2 t + \V(j)\ 2 dxdT + 2 f " f jT^n^dadT 

< / " [ <g + \V(j)\ 2 dxdT+ [ " f 4>l + \fip\ 2 dvdiodT 

Jt„-i Jr<R Jt„-! J S t 

+ f [ R(f> 2 (T,R,Uj)dLudT 

< £7[0](r„_i) + (1 + r„_!)- 1+Q (e 2 £ + Ci) . 
Now the energy inequality (21) shows that for all t < r„ 

< ^M(r) + IW» < %](t) + (1 + r)- 2+a (e 2 E + d) . 

Hence we can estimate 

(E[cf>}(T n ) - d(l + r)- 2 +«) dr < E[4>\{T n -i) + (1 + 7- n _i)- 1+Q (e 2 £ + d) , 

which, together with the fact that the sequence {r„} is dyadic, implies that 
E[4>](T n ) < r-^^Krn-i) + (1 + r„)- 2+a (e 2 £ + d) . 
In particular the energy inequality (21) shows that 

E[4>](t) <e 2 E + C 1 , Vt>0. 

Therefore 

E[4>]{Tn) < r-^MCrn-i) + (1 + t„)- 2+q (e 2 £ + d) < r" 1 (e 2 £ + d) , 
which, in turn, shows that 

E[(p]{T n ) < (1 + r„)- 2+a (e 2 £ + d) + r^ivi) 
< (1 + r„)- 2+a (e 2 E + &) + t- 2 (e 2 E + C\) 
<(l + r„)- 2 +« (e 2 E + C 1 ). 
The proposition then follows by the fact that r„ is dyadic. □ 



4 Pointwise Decay of the Solution 

A key ingredient for proving the global existence result for nonlinear problem is to derive a 
pointwise decay of the solution. Merely decay of the energy flux E[4>](t) we have obtained 
previously is not sufficient unless we show the decay of the energy for the higher derivatives 
of the solution. In Minkowski space, this is a direct consequence of the existence of global 
symmetries. On our inhomogeneous background, no such global symmetry exists. To solve this 
problem, we commute U g with £1, T and control the error terms , which are supported in the 
cylinder {(f, x)||x| < \R}, by using elliptic estimates. 

Lemma 3. The wave operator U g has the following properties: 

gVdi, =D g - g m d tt - 2g 0i d ti - -±=d a (g^V^G)d fj , (40) 
p g ,n]=f a0 d a p + fdp, (41) 



[D g ,T] = -d t g ap ■ d aP - d t ^—d a (g a W-G)) dp, (42) 

where f a ^, f 13 are smooth functions supported in the cylinder {(t,x)\\x\ < ^R} and satisfy 

\f al3 \<H, \f\<H. 
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Proof. We only have to recall the definition of the covariant wave operator 

and the fact that the metric g is a perturbation of the Minkowski metric inside the cylinder 
{r<\R}. □ 

The equalities (41) and (42) show that in order to derive the energy estimates for Tcj> or Qcj>, 
we need to estimate the error term d a pcj>. Since £lg, Tg are supported in {(t, x)\\x\ < ^R}, we 
use elliptic estimates to control the error terms. For convenience, we may omit the summation 
sign. 

Lemma 4. Suppose g satisfies the smallness condition (19) or conditions A\ and Ai- Then 
[ T2 f Id^dxdT^Eidt^ + EMT^+DidtF^l+DiF}^, (43) 

f |a Q ^| 2 ^<i?[a^](T+) + J D[9 t F]^ + J B[0](T+) + ^[F];+ fl , (44) 

Jr<R 

where t+ = max{r — R, 0} and d a fj4> * s the second order derivative in Minkowski space. 

Inequality (43) will be used when we commute U g with or T in order to obtain the energy 
decay for Q.cj> or Tcj). It shows that the error terms coming from commutations are under control. 
Inequality (44) is useful for deriving the pointwise decay of the solution when r < R. 

Proof. When g satisfies the smallness condition (19), the error term can be absorbed because of 
the smallness assumption. For the other case, we show that the inequality (43) can be reduced 
to the condition (8). 



If g satisfies condition (19), choose a smooth cut off function x{ x )> sucn that 

X{x) = 



1, |z| < \R, 
0, \x\ > R, 



and |x| < 1, |Vx| < j^, | A%| < where we denote V = (d\, 82,83). Using the inequality 

(a + b) 2 < (l+p)a 2 + (1 + I)6 2 , 

P 

we obtain 

I \d tj cp\ 2 dx< f \dM)\ 2 dx= [ \A( X <P)\ 2 dx 

Jr<±R Jr<R Jr<R 

<\ ( \A0\ 2 dx + C R [ \cb\ 2 + \V<p\ 2 dx 

< I [ I V g l] d l3 -4\ 2 dx + 6 / I V h^d t3 cj>\ 2 dx + C R [ \cj)\ 2 + \\7cj)\ 2 dx 

< \ I Ytfih.o 2 '!.,- + MH 2 [ \d l3 cj)\ 2 dx + C R f \cj)\ 2 + \\7cj)\ 2 dx, 

Z Jr<R. Jr<\R Jr<R 

where C R is a constant depending only on R. Recall that H < 700 ( 1+ °L R j a +i ■ Absorbing the 
second term, we have 

/ \d l3 cj)\ 2 dx<\ f I Y^g^d l3 cP\ 2 dx + C R f \cj>\ 2 + \Vcj>\ 2 dx. (45) 
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Then from equation (40), we have 

I 5> tf fly*l 2 < \\d a t<P\ 2 + C|D s 0| 2 + C\d<t>\ 2 . 

Add 2\d at <j>\ 2 to both sides of (45) and then integrate from n to r 2 . Then the integrated local 
energy inequality (29) implies that 

f 2 I \d a p<P\ 2 dx<4 r [ \d at cj ) \ 2 dx + C R {E[ ( j>]{T 1 )+D[FY T l). (46) 

Jn Jr<\R Jn Jr<R 

We have to bound the first term which can be estimated by the integrated local energy inequality 
for d t (f> — T(f>. Commute the equation (1) with the vector field d t . Using the identity (42), we 
obtain 

P g d t ^\ 2 < l l\Y^d t g^d a ^\ 2 + C\d t U g ^\ 2 + C\dcj>\ 2 
< 17H 2 \d aP 4>\ 2 + C^tU^ 2 + CW| 2 
for r < \R. Denote C g = 18 2 ^ 2 (1 + R) a+1 . By (30), we have 



\d at <t>\ 2 dxdT < C g E[d t 4>](n) + C g D[n g dt<f>]% 

Ti Jr<R 



< CgEidt^in) + CD[d t F] T T \ +C r f \d^\ 2 dxdT 

+ c g nH 2 (i + \-R) a+1 r [ \d a ^\ 2 dxdr 

2 Jn Jr<±R 

< CgEldtMn) + CD[d t F}ll + CE[<p}(n) + CD[F}ll 
+ C g YJH 2 {\ + ]-R) a+1 f 2 f \d ap( p\ 2 dxdT. 

1 J Tl Jr<±R 

Plug this into (46). Notice that 

o, 1 i 2 1+Q x 68 x 18 2 x (a/3) 2 
AC g nH 2 {\ + -R) a+1 < — < 1. 

We thus have shown the estimate (43) by absorbing the first term J \d a /3<j)\ 2 dxdr on the right 
hand side. 

For the case when g satisfies conditions Al and A2, Proposition 2 shows that (<jr lJ ) is uniformly 
elliptic. Thus by elliptic estimates, together with the conditon (8), we obtain 

r i \d a ^\ 2 dxdr< r [ \d at 4>\ 2 dxdT+ r [ iv^^^+a^ 

Jn Jr<\R Jn Jr<\R J n Jr<R 

< [ 2 I \d at cf>\ 2 dxdT + £7[0](n) + D[F}ll 

Jri Jr<±R 

< Eldtflin) + D[d t F]H + E[<p}( n ) + D[F]l\. 
In any case, we have proven (43). 



Now, we use (43) to prove (44). Since (g 1 -*) is uniformly elliptic as having shown in 
Proposition 2, elliptic estimates together with the estimate (40) imply that 



/ \d aP 4>\ 2 dx<( \d at 4>\ 2 dx+ [ W 3 d l3 4>\ 2 + \<b\ 2 dx 

Jr<R Jr<R J r<2R 

<E[dMr)+E[4>]{r)+ [ \F\ 2 + \d tt cf>\ 2 + \cf>\ 2 dx, (47) 

Jr<2R 
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where we recall that h al3 are supported in r < ^R. It remains to bound the integral on the 
larger ball r < 2R. We first consider the case when t > R. Take n = r — R and t 2 = t + i? in 
(20). We have 



/ / \d t <j>\ 2 + cj> 2 dxdt< / J dvol<E[ct>}(T-R) + D[F]^. 

Therefore, using Sobolev embedding, we have 

<EMt-R) + D[FY t ±%. (48) 



4> 2 dx 

r<2R 



Similarly we have 



We claim that 



\F\ 2 dx 

r<2R 



+ « ( 4 9) 



/ \du^\ 2 dx<E[dMT-R)+D[d t FY T _ R + E[4>]{T-R) + D[FY T _ R . (50) 

Jr<2R 

In fact, consider the region bounded by E t „_r and t = t. Take X = T in (9). We get 

f JlW\n^da = f JlW\n^da - f f Fd t (f) + K T [0]dvol. (51) 

Jr<2R JS T _ B n{KT) JT-RJr<R+t-T 

Notice that the metric is flat when r > \R. We can estimate 

f I \K T [<P]\dvo\< f f \d^\ 2 dxdt<E[4>](r-R) + D[F]l_ R . 

JT-RJr<R+t-T JT-RJr<R 

Using Cauchy-Schwartz inequality, we obtain 

f / \Fd t ct>\dvol< I* I |F| 2 (l + r)"+ 1 + {9 'f a+1 dxdt 

<E[4>]{t-R) + D[F]t_ r . 
In any case, g satisfies the condition A\. Hence by Proposition 2, we can show that 

/ \M\ 2 dx< f Jl^do 

Jr<2R Jr<2R 



lr<2R Jr<2R 

and 



/ J^[4>]n^da< f jT[4,]n»dv<E[4>](T-R). 



/E T _ R n{t<r} J^r-R 

Replace <j> with d t <j) in (51). Then (43) implies that 

\ \d tt <f>\ 2 dx < E[d t 4>]{r -R) + Dp g d t <PY T _ R 

Jr<2R 

<E[dt4>](T-R) + D[dtF]l_ R + [ f \d a p<f>\ 2 + \dcf>\ 2 dxdt 

Jt-R Jr<\R 

< E[d t <P}(r -R) + D[d t FY T _ R + E[cj>]{r -R) + D[F%_ r . 

The inequality (44) then follows from (47), (48), (49), (50) and (21). 

When t < R, notice that the initial data are supported in r < R. We conclude that the 
solution 4> is supported in {r < r + R}. And the above inequalities still hold when replacing 
t-R with 0. □ 
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Remark 6. The proof for the smallness assumption case shows that the condition (8) can be 
replaced by assuming that the deformation tensor tt^ is sufficiently small. This nevertheless 
still allows g to be far away from Minkowski metric. 

Having proven the above lemma, we now establish the main proposition in this section. 

Proposition 6. Suppose F satisfies the conditions in Proposition 5. Then on E r; for all a < 

5 < 1, we have 

f |r^| 2 dw<(l + r)- 1+4 (e 2 E + C 1 ), r > R, (52) 
f r\ct>\ 2 cLo < (1 + r)- 2+a (e 2 E + &) , r > R. (53) 

If in addition, 

(a) dtF satisfies the same conditions in Proposition 5; 

(b) E[d t <j>]{r) < (e 2 E + d) (1 + r)- 2 +«, 
then 

\<f>\*<(l + T)- 2+a (<?E + C 1 ), r<R. (54) 

The estimate (54) is stronger than (52) and (53) due to the extra condition (b) together 
with the robust elliptic estimates on compact region. 

Proof. Inequality (53) follows from Lemma 1 and Proposition 5. We use the p-wcighted energy 
inequalities to prove (52). Estimate (38) implies that 

/ r(d v iP) 2 dujdv< {l + T )- 1+a (e 2 E + C 1 ) , Vr > 0. (55) 
Interpolate with (34). We get 

f r 1+5 - a {d v ^) 2 dujdv< {1 + t)- 1+s (e 2 E + C 1 ) , Va < S < 1. 
Hence using (53), we can estimate 

\tp\ 2 (T, v,u>)duj < / \ip\ 2 (T,v T , cv)duj + I / / \d v ip\diodv) 

<(l + r)- 2 +« (e 2 E Q + C\)+ f V f r 1+s - a \0 v ib\ 2 dujdv f ' r'^^dv 
<(l + r)- 1+s (e 2 E + C 1 ), 
where v = and 5 > a. 



To prove (54), using(44), (15), we show that for r < R 
H 2 < / \d^\ 2 + <f> 2 dx 

Jr<R 

< E[dMr+) + D[d t FY T t R + EMt+) + D[FY+ R 
<(l + r + r 2+a (e 2 E + C 1 ) 
<(1 + t)- 2 +« (e 2 E + C 1 ), 

where t+ = max{r — R, 0}. □ 
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To obtain the poitwise decay of the solution, we need to estimate E[Q,4>](t), E[dt4>](T). 
Since n g does not commute with the vector fields 0, T, we use elliptic estimates to control the 
errors coming from commutation. We establish a proposition that gives the decay of the energy 
E[Chf>]( T ), E[d t 4>]{r). 

Proposition 7. Let the vector field X be dt or VI. Assume that F, X(F) satisfy conditions 
(a), (b) and (c) (without E[4>]{ti) on the right hand side) in Proposition 5 . 



(1) IfX = d t , then 

(2) IfX = fl and 
then 



E[dt4>]{r) < (1 + r)- 2+a (e 2 E + d) . (56) 
E[d t <P](r) < (l + r)- 2+a (e 2 E Q + C 1 ), 

E[Q<f>](T) < (1 + r)- 2+a (e 2 E + C\) . (57) 



Furthermore, in both cases U g X{4>) satisfies conditions in Proposition 5, with a new constant 
C{C\ + e 2 E ), where C does not depend on F or <p. 

Proof. Using the identities (41), (42), we can write the equation for X(<j>) 

U g X<j> = X(F) + f afi d af ^ + /%</>. 

Recall that the metric is flat when r > ^R. We infer that U g X{4>) satisfies conditions (a) and (b) 
in Proposition 5. We show that U g X(cp) also satisfies condition (c). In fact, using the estimate 
(43), we can show that 

Dp g X(<f>)]% < D[X(F)]% + r ( |<W| 2 + \d<j>\ 2 dxdr 

Jn Jr<\R 

< D[X{F)Y T \ + D\d t FY T \ + D[F}H + E[dMn) + E[<t>}(n) ' 

< (Ci + e 2 E )(l + n)- 2 ^ + £7[^](ti), 

where we have shown E[4>](t) < (C\ + e 2 Eo)(l + T )~ 2+a by Proposition 5. Therefore 

(1) when X — d t , the above inequality already implies that U g X{<f>) satisfies the condition (c). 

(2) when X = ft, the extra condition on E[d t cj)}{T) also indicates that O g X(<f>) satisfies the 
condition (c). 

Once we have verified that n g X{<j>) satisfies conditions in Proposition 5, we can conclude the 
estimates (56), (57). □ 



This Proposition shows that although D g does not commute with dt, the assumptions on 
the nonlinearity F are sufficient to obtain the decay of E[d t <f>]{T). However, to prove the decay 
of £7[Q0](r), we need to show the decay of E[d t (f>] first. The idea is that we first commute the 
equation with T. Then pass the T derivatives to fi derivatives. 

Corollary 5. Suppose Q k T 3 F satisfies the conditions in Proposition 5 for all (k,j) G A. Then 

E{Q k T^}(r) < (e 2 E + Ci) (1 + r)- 2+a . (59) 

Proof. We prove by induction on k. For k — 0, j < 8, the estimate (5) follows from the first 
case in Proposition 7. Then the second case in Proposition 7 implies that (5) holds for k = 1, 
j < 7. Repeat this again until we have arrived at the case when k — 5 and j < 3. That covers 
all the cases and the corollary follows. □ 

Remark 7. In Minkowski, we only have to commute □ with T for 3 times and with Q, for 5 
times. 
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5 Boostrap Argument 



To solve our nonlinear problem, we use the standard Picard iteration process. We prove, by a 
boostrap argument, that the nonlinear term F decays. 

Proposition 8. Suppose the nonlinearity F satisfies the following conditions 

(a) E J r<R \n k T3F\ 2 dx<2E e 2 (l+T)- 3 +«; 

(kj)eA - 

(b) E J r<R \Vn k TiF\ 2 dx<2E e 2 (l + T)- 3+a ; 

(k,j)£B - 

( c ) E It? fs \n k T3F\ 2 r 3 - a dvol<2E e 2 (l + T 1 )- 2+a . 

(k,j)eA 



Then 



I \n k T:>F\ 2 dx<E 2 6 4 (l + T)- 3+a , (60) 
(k,j)eA Jr < R 

V f \V9, k T 3 F\ 2 dx <^e 4 (l + r)- 3+Q , (61) 

„ Jr<R 



(k.j)eB 

[ 2 [ |^r J -F| 2 r 3 - Q dvol<^ e 4( 1 + ri )-2+« ; ( 62 ) 

where V denotes the covariant derivative on {t = r,r < R}. 

Estimates (60) and (62) are sufficient to conclude our theorem. The extra boostrap as- 
sumption (b) is used to prove (60). We use elliptic estimates to prove (60), which, in turn, 
implies (61). For (62), we rely on the p-weighted energy inequality and the null structure of the 
quadratic nonlinearity of F. The assumptions here, together with Proposition 5 and Proposition 
7, imply that 

E[n k T^](T)<E e 2 (l+r)- 2 + a , V(fc, j) G A. 

Since the cubic or higher order nonlinearities of F behave better, it suffices to consider the 
quadratic nonlinearity A^d^d^cj) of F, where the constants A^ v satisfy the null condition. For 
all (k,j) G A, it is known that 

Cl k T j F = ^T^A^d^dfKj)) = A af3 d a n kl T 3l <l)- dpQ k2 T j2 (f), (63) 

where notice that [T, d a ] = 0, [fi, d a ] — or dp up to a constant and ki + k 2 < k, ji + j 2 < j- For 
simplicity, in the sequel, we denote <j>\ = Q^T- 71 ^, <f> 2 — fl k2 T^ 2 (/). Before proving Proposition 8 
in details, we prove a lemma exhibiting the properties of the sets A and B. 

Lemma 5. Assume (k\ + k 2l j\ + .72) G A. Then 

(1) (h + 2,ji + 1) e A for at least one i G {1,2}; 

(2) (ki,ji) G B for at least one i G {1, 2}; 

(3) If (k,j) eA or B then (k',f) e A or B for any k' < k, f < j. 

Proof. For the first property, since k\ + ji + k 2 + j 2 < 8, without loss of generality, we assume 
k\ + ji < 4. If ki < 3, then (hi + 2, ji + 1) G A by definition; If fci > 4, then k 2 < 1 and 
k 2 + 2 + j 2 + l <4 + 3 <8. Thus (fc 2 + 2,j 2 + l) G A. 

For the second property, without loss of generality, we assume ji < j 2 . Then k\ + ji + 2 < 
k\ + ji + k 2 + j 2 < 8, which shows {k\,j\ + 2) G A. By definition, (fci,ji) G B. The third 
property holds by the definition. □ 
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5.1 Proof of (60) and (61) 

Inside the cylinder {r < i?}, it is not necessary to require the null structure of the quadratic 
nonlinearity of F. To estimate F, we rely on the robust elliptic estimates to obtain the pointwise 
bound of the solution and control the rest by the energy inequality. 
We first prove (60). Notice that 



53 \n k T>F\ 2 < \d<h\ 2 \d<h\ 2 . 

(k,j)£A 



(64) 



Since k\ + fc 2 < k, ji + 32 < j, without loss of generality, assume (k\,ji) € B according to 
Lemma 5, that is, (k\,ji + 2) e A. We claim that 



\d^\ 2 <e 2 E (l + r)- 2+a , r< l -R. 



(65) 



In fact, for dt4>\, since (fei, ji + 2) e A, Proposition 7 shows that 

E[d tt ^ 1 }<Eoe 2 (l + T)- 2 + a 
and d t O g (d t (pi) satisfies the conditions in Proposition 5. Therefore Proposition 6 implies that 

\d t <j>i\ <e 2 E (l + r)- 2+a , r<±R. 
For |V0i| 2 , using elliptic estimates and the inequality (44), we can show that 

\vd aP <h\ 2 dx < f |^'%v^i| 2 + iva^^i 2 + id^i 2 + \d<f>\ 2 dx 

r<ifl Jr<R 



< 



< 



J (□„ - g 00 d tt - 2g 0l d u - -jL=d a (g a PV=G)d^ V0i 

+ / \vdpdt(h\ 2 + \d a p(i>i\ 2 + \d(f>\ 2 dx 

Jr<R 

f \Vn g 4 1 \ 2 + \Vd tt <j>i\ 2 + \Vcb 1 \ 2 dx+ f Id^faf + ldfafdx 



lr<R 
+ 



< 



f \vdpdtfal 2 + \d a ^ 1 \ 2 dx + \d<f>\ 2 

f \Va g( j) 1 \ 2 dx + 6 2 E (l + T+)- 2 + a , 
where t + = max{r — R, 0}. Now can write 

k<ki,j<ji 
k+j<ki+ji 

+ E F^d aP n k T^ + F^d a n k T^, 

k<ki,j<ji_ 

where fZf,f%,F?f are smooth functions depending on g and are supported in r < \R. 
Therefore, using the boostrap assumptions (b) and the estimate (44), we obtain 

/ \Vd a ^ 1 \ 2 dx<Eoe 2 (l+T)- 2+a + V f \V8 a ^ v \ 2 dx, 

where V < 1 means ky < k\, jy < j\ and ky + jy < fci + ji- By induction, we get 

/ \Vd afj cj )1 \ 2 dx<E e 2 (l + T)- 2 + a . 

Jr<hR 
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Then Sobolev embedding and the estimate (44) imply that 



! W Jr<h 



C5 ( B 1 K ) Jr<kR 



< I WiifaP + ldaBfaf + Wifdx (66) 



2C1 _l_^-2+a 



<E e'(l + T)- 

We thus have shown the desired estimate (65), which implies that 



V f \n k T^F\ 2 dx< (l + T)- 2+a E e 2 [ \dct> 2 \ 2 dx 

<(l + T)- 2+a E e 2 E[fo](T) 
< (l + T)- 3+a E 2 e 4 . 

When 7}R < r < R, we use the angular momentum f2. Observe that k\ + k 2 < 5. Using Sobolev 
embedding on the unit sphere, we always have 

/ \d<h\ 2 ■ \d<h\ 2 du < [ \d<j>v\*du- I \d<t>2>\ 2 du, (67) 

where ky < k\ + 2, k 2 ' = k 2 if k\ < k 2 - otherwise, ky = k\, k 2 > < k 2 + 2. In any case, 
(ky,j\), (k 2 ',j 2 ) G A. Since j\ + j 2 = j, without loss of generality, we assume (ky,ji + 1) G A. 
Therefore, we have 



(68) 



f \dcj)y\ 2 duj< [ \d(j)y\ 2 + \\7d<f)y\ 2 dx 
Jul Jr<R 

< E[<f>y](r) + E[d t M{r + ) + D[d t FY+ R + E[cj>y]{r + ) + D[FY+ R 
<(l + T)- 2 + a E Q e 2 , 

where t + = max{r — R, 0}. Then (67) implies that 

V I \n k T^F\ 2 dx< (l + T)- 2+a E e 2 [ \8cj) 2 ,\ 2 dx 

(k,j)eA J i R < r < R Jr < R 

<{l + T)- 2+a E e 2 E[^]{r) 

<(l + r)-^Ele\ 

This completes the proof for (60). 

Remark 8. We remark here that (68) is only true when r is bigger than a constant. That is 
why we need to distinguish the two cases: r<\R and\R<r < R. 

Now we use (60) to prove (61). By (63) and (67), for all (k,j) G B, we can show that 

f \vn k r^F\ 2 duj< f \d^\ 2 - \d^<f> 2 \ 2 du < [ \d<t>y\ 2 dij- [ \d M <t> 2l \ 2 dLo, 

where ky < k\ + 2, k 2 > = k 2 if k\ < k 2 - otherwise, ky = k\, k 2 > < k 2 + 2. In any case, 
(ky,ji), (k 2 ',j 2 ) G B. Thus by (65) and (68), we have 

j \dM 2 du<(l + r)- 2+a e 2 E , r < R. 

On the other hand, (44) implies that 

\d Pl cj )2 ,\ 2 dx<(l + T)- 2 + a 6 2 Eo. 



I I 

Jr<R 
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Therefore 

-3+Q 



V / \VSl k T*F\ 2 dx < (1 + T)- 2+a e 2 E f \ 2 dx < E 2 e 4 (l + r)" 

(k,MB Jr < R Jr < R 

We thus have shown (61). 

Remark 9. It is not necessary to require the nonlinearity to satisfy the null condition when 
r < R. Moreover, when r < R, the nonlinearity can be any form of 

F = F(d<f>,4>), F(0,0) = 0, DF(0,0) = 0. 

5.2 Proof of (62) 

For this part, we rely on the p-weighted energy inequality and the null structure of the nonlin- 
earity. For the null form N = A^d^d^cf), we have 

Q k TiN{d<f>, d(f>) = 

where (pi — Q, kl T^(j), <j> 2 — fl k2 T^ 2 (f) and k\ + k 2 < k,ji + j 2 < j. The p-weighted energy 
inequality is an estimate in terms of ip — r<f> instead of p. For this reason, we expand F in ip. 

Lemma 6. If N is a null form, then 

r 4 \n k T>N\ 2 < ]T <t>\<t>\ + 0? ' r 2 d 2 r <p 2 + \fH 2 \fH 2 + \d v H 2 \duH 2 , (69) 

1,2 

where v = , u = . 
Proof. In fact, notice that 

r 2 N{d<f> u d<h) = 0102 + r(<^ 2 ) r + N{d^d^ 2 ) 

and 

\N{d^,d^ 2 )\ < \a v ih\ ■ \duH + \d u H ■ \d v ip 2 \ + |yvi| • \?M- 

The lemma then follows. □ 

Using this lemma, it suffices to consider the four terms on the right hand side of (69). We 
handle the first three terms in a uniform way. Let $i be <p\ or fipi; $2 be <p 2 , rd r (j) 2 and y-ip 2 
correspondingly. Since k\ + k 2 < 5, using Sobolev embedding on the unit sphere as we did in 
(67), we obtain 

f |$i| 2 |$ 2 | 2 ^< f |$i,| 2 dw- f \$ 2 ,\ 2 du), (70) 

where k v < k\ + 2, k 2 > = k 2 if k\ < k 2 ; otherwise, ky = ki, k 2 i < k 2 + 2. Here we omit the 
summation sign. For the third case when $1 = yVii ^2 = ^ip2, without loss of generality, we 
assume k\ < k 2 . Thus by the fact that y/ipi = Q<p\, we can always write $1/ as Q k T : >cj> for some 
(k, j) € A. Then (52) shows that 



r 2 



f \<P v \ 2 du <(l + T)- 1+s e 2 E , V5>a. 

J UJ 



Therefore, by the integrated energy inequality (20), we have 

(•T2 P PT2 poo 



pT-2 P PT2 POO P 

/ / r- 1 - a $ 2 1 <S> 2 2 dvo\= / / / r 1 - a <S> 2 1 <$> 2 2 dvdojdT 

J n J St J Tl J V t J UJ 

pT2 POO P P 

< / / r 1 -* / \<$> v \ 2 duj / \<$> 2 ,\ 2 dujdvdT 

J T\ J V T J OJ J UJ 

<(l + n)-^r/ Mdvol 

J Tl J St 



(71) 



< (1 + T 1 )- 1+5 e 2 E (l + T 1 )- 2+a e 2 E 

<(l + n)- 2 +«eX. 
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It remains to estimate the main terms \d v i/j\\ 2 \d u ip 2 \ 2 in (69). There are two cases according 
to which one is bigger: k\ or fc 2 . 

When ki < k2 1 in particular we have ki < 2. The idea is that we bound |9 t ,V'i| uniformly 
and then control \d u 4>2\ 2 by the energy flux through the null hypersurface v = constant. Wc 
first establish a lemma to show that the energy flux through v = constant is bounded. 

Lemma 7. Consider the region D = [tt 1,1x2] x [vi,oo) cS,x [ti,t 2 ]. Then 

/ / {du^fdwdu < (1 + Tl )- 2+a e 2 E a . 

J Ui J \oj\ = 1 

Proof. Back to the energy equation (9), take X = T on the region D. We get 

J 2 jT[<j> 2 ]n»do-+ J jT[fo}n^da = J J? [foK da + jj[£j]n"dt7 



V>V\ V>V\ 
U—Ui U—U2 



^g4>2 ■ d t 4>2dvo\- 
D 



Thus applying Cauchy-Schwartz inequality to the last term, we obtain 
/ r 2 (d u (j)2) 2 dLjdu<2 J^[<f>2}n"da 

< E[<h](n) + D[n g <h]Z + ^ J^Hll_dvol 

<(f + Ti)- 2 +«e 2 i?o. 

Then (53) implies that 

I j (d u ip 2 ) 2 duidu = / / r 2 (d u (f>2) 2 + d u (r(f>l)dujdu 
= / r 2 (d u 4>2) 2 dudu + / rcj)\duj 

J u\ J — 1 J to Ul 

< (1 + T 1 )- 2+a e 2 E a + (1 + T 2 )- 2+a e 2 E 
<(l + T 1 )- 2+a e 2 E . 

We continue our proof when k\ < 2. The above lemma shows that 

poo pu(v) 



□ 



/ 2 / r-^I^ViH^pdvol^ / / |a ti ^ 2 | 2 r 1 - a |a i ,Vi| 2 d«dwdw 

rOO pu(v) 

< / \d u ifa\ 2 dwdu-axipr 1 - a \d v il>i\ 2 dv (72) 

Jv T1 Ju T1 u, u 



POO 

<(l + T 1 )- 2+a e 2 E / supr^l^Vil 2 ^- 

Jv T1 »,« 

We use Sobolev embedding to estimate supr 1_Q |9„'i/'i| 2 . First, on the unit sphere, we have 



I \^ a d v ^\ 2 dLo <J2[ \d v n a in\ 2 u = f \d v ip v \ 2 duj, (73) 



I 2 < 

a<2" UJ a<2 
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where ky < k\ + 2 and we omit the summation sign. Then on [m, u 2 ], we have 



+ / r 1 - a (d u d v ^ 1 ,fdu+ / r- a (d v M 2 du 

+ / ^'"(^ipy) 2 + r 3 - a \Fy\ 2 du, 

where we use the wave equation (32) and U\ = u T = ^^,U2 = u(v). Integrate over the unit 
sphere. We can show that 

/ / r^ldvipvfdcudv < f r 1 - a {d v ^y) 2 dvdoo 

J v T1 J uj J St 1 

+ I 2 I ^-"(d^v) 2 + {yVttPyfr 1 -* +r 3 - a \F 1 ,\ 2 dvdiudT 1 

where f = ^. We claim that the above inequality can be bounded by a multiple of e 2 E . 
In fact, the first term can be bounded by (1 + ti)~ 1+q e 2 Eq by (55); the second term can be 
bounded by e 2 E by (35); the third term can be controlled by (l + Ti)~ 2+a e 2 E by the integrated 
energy inequality (20) (notice that ky < k\ + 2 < 4); The last term is good by our bootstrap 
assumptions and Proposition 7. Summarizing, we have shown 



f 



supr^ldyipxfdv < e 2 E . 

Jv T1 u,u> 

Plug this into (72). We get 

P r 1 "" / \d v ^\ 2 \d u ^ 2 \ 2 dvdojdT < (1 + T X y 2 + a ^E 2 . 

When k2 < fci, that is k 2 < 2. For this case, we control Ic^il 2 by the p-weighted energy 
inequality and bound d u ip2 uniformly. First using Sobolev embedding on the unit sphere, we 
have 

\duH 2 < E / 1^ ^ duj ~ E / l$Wi| 2 dw - / \du^\ 2 du, (74) 
where ky < k 2 + 2 and we omit the summation sign. Therefore, by (35), we obtain 



T2 



r 



l-a 



d v fpi | 2 \d u ip2 \ 2 dvcLodT 



n JS 

t>T 2 



< ( * f r^ld^l 2 ■ f \d u i, T \ 2 duj dvdwdT (75) 
<(l + r 1 )- 1+a e 2 E f T \upr- a I \d u ^\ 2 dudT. 

Jti v Jul 
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Then for all v G [v T ,oo), we can show that 

r~ a (du^) 2 < r- a (d u ip 2 ,) 2 \ + 



r^- a \d u ^\ z dv 



+ 2 



/ r~ a \d u ip 2 > ■ d v d u tp 2 '\dv 

J V 



/>oo 

< r- a (d u iP 2 ,) 2 \ v=VT + / r-^du^dv 

r oo /*oo 

+ / r- 1 - a (d u ij 2 ,) 2 dv+ / ^{dvdu^.fdv 

< r (d u ^) \ v=Vt2 + ^ rl+a dv 



+ / r 1 - a (^ 2 ,) 2 dv+ / r^ a \F 2 ,\ 2 dv. 
Intcgrating over the unit sphere, we conclude that 

[\upr- a f \d u i> 2 >\ 2 dLodT< f 2 f r- a (d u ^ 2 ,) 2 \ v=VT dr 

+ T l {du f+f +r 1 - a (W 2 ,) 2 +r 3 - a \F 2 ,\ 2 dvdLudT. (76) 

Jr, JSt r 



In JS. 

Reparametrize the first term. We get 

^ f r- a {d u i> 2 ,f\ v=VT dr= f ^ f r- a \d u 4> 2 ,\ 2 dLodu< f ^ f \d u ^\ 2 dujdu, 

which can be bounded by (1 + Ti)~ 2+a e 2 E by Lemma 7. The second and third term in (76) 
can also be bounded by (1 + Ti)~ 2+a e 2 E because of (20) (notice that k 2 > < k 2 + 2 < 4). The 
last term in (76) is good by our assumptions and Proposition (7). Therefore 

f 2 I r^d^^d^dvdujdr < (l + r 1 )- 1+a e 2 E ■ (1 + r 1 )- 2+a e 2 E < (1 + Tl )- 2+a e 4 E 2 . 

This concludes the proof of (62) and hence Proposition 8 follows. 

6 Proof of the Main Theorem 

We used the foliation £ T , part of which is null, in the previous argument. However, we do not 
have a local existence result with respect to the foliation E T . We thus use the standard Picard 
iteration process. Take 4>~i(t,x) = 0. We solve the following linear wave equation recursively 

I U g(t,x)4>n+1 = F((j) n ,d(j) n ), ^ 

\<f) n+ i(0,x) = e<j) O (x),d t (j) n+1 (0,x) = efaix). 

Now suppose the implicit constant in Proposition 8 is C\, which, according to our notation, 
depends only on R, a, A, h and Co- Set 

1 



V C\Eq 

Then for all e < e , we have 

C\(^E 2 < € 2 Eq. 
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Thus by the continuity of F(<j> n ,d(j) n ), we in fact have shown that the nonlinear term F satisfies 

I \n k TjF{<j> n ,dcj) n )\ 2 dx<C 1 EZe 4 {l + T)- 3+a <E e 2 (l + r)- 3+a , Vn, 
( k ,j)eA Jr < R 

f 2 I |r> fe T^F(0„,a ( /»„)|V- Q dvol<S o e 2 (l + T)- 2 + Q , Vn. 

Then Corollary 5 implies that 

E[n k T^ n ](r) < E e 2 (l + T)- 2+a , V(fc,j) € A. 

We remark here that all the implicit constants are independent of n. 

Then Proposition 6, together with the Sobolev embedding on the unit sphere, implies that 



^ |ft fe T^„| < v^e(l + r)-'(l + |t-r + i2|)- 1+ ' a , 

k<2,j<2 

\K\ <s \fE~ e{l + r)- 1 ^ + |t - r + R\)~^ s , Vtf > a. 

We still need to show that the solution is C 2 . The first step is to show that <j> n is uniformly 
bounded in C 1 . When r < the inequality (65) implies that 

\dV, k T 3 4> n \ < V^e(l + r)-5(l + T )- 1+ 5 a , Vfc < l,j < 2. 

When ii? < r < i?, we use (68) and Sobolev embedding on the unit sphere and we can obtain 
the same estimates. When r > R, noticing that ^ = ^, it suffices to consider d r (j) n . First the 
inequality (76) shows that 

f T2 snpr- a f \d u M 2 + \dtd u M 2 dujdT<(l + T 1 )- 2+a e 2 E^ 

Jti v Juj 

where ip = ril k Ti(j> n , k < 3, j < 2. Using Sobolev embedding on S 2 x [ti,T2], we obtain 

\rd u n k T^ n \ 2 < <f> 2 +r a (l + T)- 2+a e 2 E , Vfc < l,j < 2. 

Since d u = d t - d r and |O fc T^„| 2 , \d t 9, k T^ n \ 2 < (1 + r)" 1 ^ + r)- 2 + Q e 2 £;o and a < 1, we 
obtain 

\d r n k T^ n \ < (1 + r)"3(l + r)- 1+ ^/^e, Vfc < 1, j < 2. 
In particular, we have shown that 

\3n k T^ n \ < (1 + r)-*(l + T)- 1+ ^v^e, Vfc < l,j < 2. 

To show that <fi n is also bounded in C 2 . Outside the cylinder {r < jR}, we use the equation. 
Back to the equation (77), we can represent d rr (j) n as follows 

g rr d rr <t> n = F^n-i, dfyn-x) - g afi d a0 cj) n + g rr d rr cj) n - y=^d a (g^^G) 

Since we have shown that 

|5M |tt 2 0„|, \d tt K\, IM Wu\, W4> n \ < (1 + r)"'(l + ri)- 1+ 3 Q V^C Vn, 
we thus can estimate 



\d 2 K\< Y. m k T^ n \ + \d rr 4n\<^E~oe(l + r)-Hl + T)- 1+k i< 

k,j<i 
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Inside the cylinder {r < jR}, we rely on elliptic theory. First, we have the elliptic equation for 



g ij d ij( t> n+1 = F(<f> n ,d<p n ) - ~^d a (g af} ^G) d^ n+1 - 5 00 T 2 ^+i - 2g 0i diT<f> n+1 . 

As we have shown from estimates (66) that 

\\dn k T^ n \\ ci < \\Wn k T^ n \\ ci +||O fe T^Vn||cc < V^e(l + r)"'(l + r)- 1+ ^ 

h R i R 

for all k < l,j < 2, we conclude that the right hand side of the above elliptic equation is 
uniformly bounded in C^(Bi R ). Hence Schauder estimates [?] imply that 

ll^+ill^W <^(l + r)-5(l + r)- 1+ ^. 

A R 



In particular, the above argument shows that 



E 1^1 £ V^c(l + r)-'(l + r)- 1+ ^. 

I/3|<2 



(78) 



Now the classical local theory shows that there exists a time t* > and a unique smooth solution 
(f)(t 7 x) <G C°°([0, t*) x K 3 ) of the equation (1). Moreover, the proof of the local theory indicates 
that 

<p n (t,x) ->• (f>(t,x) 

in C°°([0,t*) x M 3 ) since the proof of the local existence result relies on the Picard iteration. 
Therefore by (78), we have pointwise bound for the solution <p 

< /^e(l + r)-5(l + r)- 1+ ^, V(t, a;) G [0, T) x R 3 . 

|/3|<2 

By a theorem of Hormander [7] that as long as the solution is bounded up to the second 
order derivatives, the solution exists globally. That is there exists a unique global solution 
(j){t,x) <G C°°(M 3+1 ) which solves (1). Moreover 



<f) n (t, x) -> <j)(t, x), (t, x) e R 3+1 . 
Therefore (f> admits all the estimates of <j) n obtained above . 
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